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Abstract 

Within the framework of the improved quantum molecular dynamics model, the medium mod¬ 
ifications on the free nucleon-nucleon elastic cross sections are investigated. By using various 
in-medium nucleon-nucleon elastic cross sections in the model, the nucleon-induced reactions on 
various targets are simulated, and the excitation functions of reaction cross sections in the energy 
range from 25 MeV to 1 GeV are calculated. By comparing the calculations with the experi¬ 
mental data, an isospin, density, and momentum-dependence medium correction factor on free 
nucleon-nucleon elastic cross sections is determined. 
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I. INTRODUCTION 


With more and more upcoming experimental data of rare isotope heavy ion collision at 
intermediate energies, people have opportunity to explore the properties of nuclear matter 
with high density and large isospin asymmetry by the tools of transport model. Both the 
mean held and the in-medium two-body scattering cross sections, coming from the same 
nucleon-nucleon (NN) interaction [l-4|, should be treated self-consistently in the transport 
model. But it is very difficult to solve the dynamical equation and the G-matrix simultane¬ 
ously, the most of transport models deal with the mean held and collision separately. So the 
accurate in-medium nucleon-nucleon cross sections (NNCS) are required by transport model. 
On the other hand, the in-medium NNCS can provide information which is very important 
to study the structure of nuclei especially rare isotopes. Also the in-medium NNCS are of 
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interest for their own sake, as they are underly related to the viscositw 
of nucleon in nuclear matter, and other nuclear transport coefficients 5 

Unlike the NNCS in the free space which can be directly measured by experiment, the 
in-medium NNCS can only be calculated by theory, such as the Brueckner-Hartree-Fock 
theory, the Dirac-Brueckner-Hartree-Fock theory, the Green function approach, and the 
relativistic mean held model. The in-medium NNCS depend on several quantities, such as 
density and isospin-asymmetry of nuclear matter, total momentum of the nucleon pair in the 
nuclear matter rest frame, the relative momentum of the nucleon pair, even the temperature. 
There have been many studies demonstrate that the NNCS in medium should be suppressed 


compared to free NNCS 2l. l8Ml9|. However, the modihcation factors given from these works 


differ signihcantly, due to the different treatment on the dependence with different levels of 
approximations. 

The other way to get the information of the in-medium NNCS is to compare the simulation 
results from transport theory with experimental data of nuclear reactions. As the basic 


NNCS are widely used. Such as the parameterizations given by G. Li 


20N22|. the phenomenological in-medium NNCS scaled by effective mass 
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where fi and p* are the reduced masses of colliding nucleon pairs in the vacuum and medium. 
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respectively, and the empirical relation 


= (1 + hP/Po)^^NN^ 


( 2 ) 


where r] is an adjnstable parameter. Consistent with theor y re snlts, the most of experimental 


evidences, snch as balance energy [27H29|, stopping power |30|, collective flow 
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etc, snpport a rednced in-medinm nncleon-nncleon elastic cross sections (NNECS). Bnt the 
reduced factor has not been determined very dehnitely. Since both the mean held and the 
two body collisions are convolnted, model dependent treatment of the nnclear potential as 
well as the collisions may yield considerable discrepancies in the model ontpnts even for the 
simplest box calcnlations. Compared to the heavy ion collision, some kind of direct reactions 
involve less degrees of freedom in the reaction process may rednce the model dependence. 
So the nncleon-indnced reaction is more snitable to stndy the in-medinm NNCS, not only 
becanse of its simple mechanism bnt also becanse it rehects jnst the in-medinm NNCS 


aronnd the incident beam energy and below the satnration density 32 h36|. 


In this work, with the improved qnantnm molecnlar dynamics (ImQMD) model, we try to 
extract the in-medinm NNECS from the experimental data of the nncleon-indnced reaction 
cross sections (RCSs). The paper is organized as follows. In Sec. II, we briehy introdnce the 
model we adopted. In Sec. Ill, we introdnce the form of medinm modihcation on NNCS, 
and the way to determine the modihcation factor. Finally a brief snmmary is given in Sec. 
IV. 


II. MODEL 


Within the ImQMDOb model 30,1371], the mean held acting on nncleon wave packets are 
derived from an energy fnnctional where the potential energy U inclndes the fnll Skyrme 
potential energy with jnst the spin-orbit term omitted 


U U\qq -|- Ccoui) 


here, Vcoui is the Conlomb energy, and U\oc = / hioc[p(^’)]dr. V\oc[p{r)] reads 
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where p, pn, Pp are the nucleon, neutron, and proton density, 6 = {pn — Pp)/{pn + Pn) i s the 
isospin asymmetry. The coefficients utilized in equation (jl]) can be transcribed, see 30|, 1^ , 
onto those usually specihed for the Skyrme interactions. 

The Pauli principle is considered carefully in the collision term. The neutrons and protons 
are treated as different particles. According to the uncertainty principle, the final states of 
two particles are required to satisfy the relation 
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The Tij and pij are the distances between the centers of wave packets of two nucleons in 
coordinate and momentum space. The possibility being blocked at the hnal state i and j 
for each of the two scattering nucleons is calculated by 


-Pblock — 1 — (1 — -Pi) (1 — Pj)- 


( 6 ) 


The Pi and Pi are the probabilities for state i and j being occupied, which can be calculated 
by 

A 
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III. RESULTS AND DISCUSSION 


The RCS for nucleon-induced reactions is calculated by the formula 

rbmax 

aR= 27rPinei(fe)M&, (8) 

Jo 

where Pinei(^) represents the probability of the inelastic scattering event at the impact pa¬ 
rameter b. 6max is the largest impact parameter at which there is no more inelastic scattering 
event taking place. Experimentally, ctr is indirectly measured by measuring the elastic scat¬ 
tering cross sections. The elastic scattering event is dehned as event that the emitted nucleon 
has its energy close to the incident energy, without considering how many NN collisions are 
happening in the event. According to the analyzing method used in experiments, here 
Pinei(^) is calculated by Pjnei = 1 — Pun — -Pel, where Pun is the probability of the events 
that the incident nucleon only passes through without changing its momentum direction 
and magnitude, namely, the unaffected events and Pei is the probability of the events that 
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the incident proton only changes its momentum direction but not its momentum magnitude, 
namely the elastic events, respectively. 

As mentioned above, the mean held and NN collision affect the nuclear reaction process 
together. The effect of the mean field on RCS should be checked firstly. As a test, four 
Skyrme parameter sets MSLO 381, BSkl 39], SkP ^|, SkI2 41] are adopted to simulate 
the reactions of p+^^Fe. In the calculations, the free isospin dependent NNCS provided by 


Cugnon et al. 


42j is adopted. The comparison between the calculated excitation functions of 


RCS and the experimental data is illustrated in Fig. [T] One can see that, with various EOS 
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FIG. 1: (Color online) Excitation functions of RCS for p+®®Fe calculated with free NNECS and 


four Skyrme parameter sets compared with experimental data. Data are taken form Ref. 


43|. 


adopted, the calculated RCSs are different. With energy increasing the difference becomes 
smaller, due to the NN collision dominates the reaction. For the case of energies higher 
than 500 MeV, the results calculated with free NNECS can reproduce the experimental 
data well. But all calculation results overestimate the experimental data below 500 MeV. 
It means that the in-medium NNECS below 500 MeV should be smaller than free NNECS 
and the in-medium NNECS above 500 MeV should be very close to free NNECS. To obtain 
correct medium correction on NNECS, the reasonable EOS tested by other method should 
be used in simulations. MSLO, one of Skyrme parameter sets which best satisfy the current 
understanding of the physics of nuclear matter over a wide range of applications 4^ , is used 
in the following calculations. 

According to the above analysis and the other studies on NNSC, in this work, the free 
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NNCS is modified according to 


,free 


+ a:, = aT + F{u,6,p)aT- 


(9) 


Where the Ugi®® and are the free isospin dependent elastic and inelastic cross sections, 
respectively. The form of F{u,6,p) is as the same as that proposed by Q. Li et al. in 


Refs. 


20 


22l |. The medium correction factor F = F^ ■ F^ depends on the nuclear-reduced 


density u = p/po, the isospin-asymmetry 6 = {pn — Pp)/{pn + Pp) and the momentum. Where 

= A-1-(1 - A) exp(-M/C), (10) 


Fs = I - TijA{u)6, A{u) = 


0.85 


( 11 ) 


1 + 3.25m 

When i = j = n, Tij = —1; i = j = p, Tij = +1] i ^ j, Tij = 0. The Fj and F^ factors are 
expressed in one formula, 


Flu = 


fo, 


F's,u—fo 


Pnn > 1 GeV/c, 
Pnn < 1 GeV/c, 


( 12 ) 


1+(pnn/po)'‘ 

with pnn being the relative momentum in the NN center-of-mass system. By varying the 
parameters A, (, /o, po and k one can obtain different medium correction on NNEGS. The 
parameter sets used in this work are listed in Table [I] and Ull Among these parameter sets, 
FUl-3 and FPl-5 are taken from Refs. 


20 


22| , FU4 and FP6 are obtained by this work. The 


TABLE I: Parameter sets used for the density-dependent correction factor F^. 


Set 

A 

C 

FUl 

1/3 

0.54568 

FU2 

1/4 

0.54568 

FU3 

1/6 

1/3 

FU4 

1/5 

0.45 


medium correction factor Fu as function of reduced density u and momentum dependence 
of Fu with FU4 at p/po = 0.5 are represented in Fig. El From the figure, one can see the 
rough character of NNFGS, with density increasing in-medium NNFGS decrease. And for a 
certain set of density dependence, with momentum increasing the in-medium NNFGS might 
be enhanced. 
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TABLE II: Parameter sets used for the momentum dependence of correction factor F^^s- 


Set 

fo 

Po [GeV/c] 

K 

FPl 

1 

0.425 

5 

FP2 

1 

0.225 

3 

FP3 

1 

0.625 

8 

FP4 

1 

0.3 

8 

FP5 

1 

0.34 

12 

FP6 

1 

0.725 

10 



FIG. 2: (Color online) (a) The density dependence of and (b) the momentum dependence of 
Fu with FU4 at p/po = 0.5, respectively. 

The various in-medium NNECS obtained from combinations by parameterizations FUl, 
FU2, FU3 and FPl, FP2, FP3 are tested by the excitation function of RCS for p-|-®®Fe. By 
using the in-medium NNECS, the descriptions on the excitation function of RCS are great 
improved. Especially the combinations of FU2-I-FP3 and FU3-I-FP3 give the best two results 
except a little deviation at low energies. According to the momentum dependence of the 
in-medium NNECS shown in Fig. [2](b), the enhancement effect of momentum correction 
dose not yet appear at such low energies if FP3 is adopted. As a test, the excitation 
function of RCS calculated by FU3 and without momentum correction is shown in subhgure 
(c) marked by FU3FP0. One can see that the momentum correction dose not obviously 
affects the results below 50 MeV, while the RCSs at high energies are underestimated if the 
momentum correction is absent. It proves our assumption and means the depressive effect 
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FIG. 3: (Color online) Excitation functions of RCS for p+®®Fe calculated with free NNECS and 


various in-medium NNECS compared with experimental data, respectively. Experimental data are 


taken form Ref. 


43| 


provided by FU3 is a little strong while the one provided by FU2 is a little weak. So we 
propose a parameter set, namely FU4, to provides a reasonable correction effect between 
ones given by FU2 and FU3. For the density correction parameter set FU4, corresponding 
momentum correction parameter FP6 is obtained by fitting the experimental data of RCSs 
at higher energies. 

More experimental data are used to test the obtained in-medium NNECS FU4FP6. From 
light to heavy targets, reactions of proton-induced on ^^Al, ®°Zr, ^^®Sn, ^°®Pb 

are simulated. The excitation functions of RCS for p-|-A calculated with the in-medium 
NNECS compared with the experimental data are presented in Fig. 01 One can see that, 
all experimental data can be quite well reproduced. It demonstrates that the in-medium 
NNECS given by FU4FP6 is reasonable. Especially, not only the excitation functions of 
RCS for the targets along the /9-stable line can be well described, the excitation function 
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FIG. 4: (Color online) Excitation functions of RCS for p+\ calculated with FU4FP6 compared 


with experimental data, respectively. Results calculated with free NNECS are also represented for 

45-48| for C, 45] for Ca, 3, 46-48] 


reference. Experimental data are taken form Refs. 


for A1 




n 


for Sn, [43|, 1481] for Zr, and (43|, l45|, 


43 , 


46l | for Pb targets, respectively. 


of RCS for ^®Ca, which is far away from the /9-stable line, can also be well reprodnced. 
Beside the proton-indnced reactions, the nentron-indnced reactions are also used to test the 
in-medium NNECS. As the results presented in Fig. [5l the excitation functions of RCS for 



EIG. 5: (Color online 
taken form Refs. 


Same as Eig. 


491-15J for C, 


51- 


but for neutron-induced reactions. Experimental data are 

d, 


55l | for Cu, and 


51 M 55 I] for Pb, respectively. 


n-l-A calculated with FU4FP6 are in good agreement with the experimental data. It means 
that the medium correction on isospin dependence of NNECS is considered reasonably. 


9 


























































IV. CONCLUSION 


By using the improved quantum molecular dynamics model, the nucleon-induced reac¬ 
tions on various targets with incident energies from 25 MeV to 1 GeV are investigated. The 
reaction cross section is found to be very sensitive to the nucleon-nucleon cross sections. By 
comparing the excitation function of reaction cross sections between ImQMD model calcula¬ 
tions and experimental data, the medium modihcations on the free nucleon-nucleon elastic 
cross sections are investigated. An isospin, density, and momentum-dependence medium 
correction factor on free nucleon-nucleon elastic cross sections is obtained. The parame¬ 
terized formula is simple to use and gives reliable results for the nuclear reactions systems 
below saturation density. 
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